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Abstract 

The dimension-four genuine anomalous quartic couplings are studied in processes 
of six-fermion production via e + e~ collisions. Complete tree-level electroweak cal- 
culations are performed including initial-state-radiation and beamstrahlung. The 
analysis of final-state distributions can be used to find kinematical cuts to enhance 
the effects of anomalous couplings. For the parameters of the custodial-symmetry- 
conserving anomalous couplings a sensitivity in the range between 1CP 3 and 1CP 2 
can be expected at 1 TeV. 
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1 Introduction 



The experimental measurement of the gauge self-couplings is an important test 
of the Standard Model (SM), that is still at the beginning at present colliders. 
The trilinear gauge couplings have recently become a subject of studies at LEP 
II and Tevatron (for recent results, see [1]), and will also be measured at the 
Linear Collider (LC) [2] . The quartic couplings are very weakly constrained by 
present experimental data, through loop diagrams [3], while the observation 
of quartic gauge coupling effects at tree-level at LEP is rather difficult when 
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photons are involved [4] and is completely outside the reach of LEP when only 
massive gauge bosons are present. 

The determination of quadrilinear gauge couplings is of particular interest in 
connection with the problem of the electroweak symmetry breaking. Indeed, if 
a Higgs boson lighter than 1 TeV does not exist, the tree-level amplitudes for 
the gauge boson scattering in the Standard Model are known to violate the 
unitarity limits at tree-level, indicating the existence of strong interactions. 
This scenario can be the result of different models of symmetry breaking and 
may be described in the general framework of effective lagrangians, where the 
electroweak interactions in the low-energy limit are represented by one general 
parameterization [5] . Different values of the parameters correspond to different 
models. The effective lagrangian is organized as a power series in p/inv, where 
p is the scale of momenta of the phenomena under study and v ~ 246 GeV 
is the scale of symmetry breaking. The lowest-order terms of this expansion 
are model- independent. The higher order terms involve new gauge coupling 
structures that give rise to anomalous vertices and contain coefficients that 
are model-dependent. 

It is worth noticing that the lowest-order non-trivial quartic couplings in the 
effective lagrangian (that are the 0(p 4 ) or dimension-four terms) involve only 
massive gauge bosons and thus their effects can be observed either in loop 
contributions, or at tree-level in processes with at least six fermions in the 
final state. 

The measurement of the parameters describing the gauge self-couplings, that 
gives the opportunity of constraining some possible models of new physics, 
will be an important objective of the LC, where significant improvements will 
be possible with respect to present colliders for various reasons. In the first 
place, the cm. energy is sufficiently high for the production of up to three 
gauge bosons in the final state, making it possible to study the quadrilinear 
couplings at tree-level. Moreover, since the deviations from the SM values of 
these couplings destroy the unitarity cancellations, their effect is enhanced at 
high energy. Finally, the special features of the LC, such as the high luminosity 
and the possibility of having polarized e + and e~ beams will provide a very 
high sensitivity to the anomalous couplings. 

Several phenomenological studies have been performed on the possibility of 
constraining the anomalous quartic couplings at e + e _ colliders in the approxi- 
mation of real vector bosons in the final state, both for the couplings including 
photons [4] and for those involving gauge bosons only [6,7]. In the latter case, 
however, for more realistic predictions, as observed above, processes with at 
least six fermions in the final state have to be examined. This kind of processes 
also allows for the analysis of final-state distributions that are not accessible 
in the real approximation, and that can give useful information. 
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The objective of the present study is to analyse a class of six-fermion (6/) 
processes where the genuine quartic gauge couplings of dimension four are 
involved. The results that will be shown have been obtained by means of 
complete tree-level calculations in the framework of the chiral approach to 
electroweak interactions. 

In Section 2 the theoretical framework and some technical details of the cal- 
culation are explained. In Section 3 the numerical results are presented and 
discussed, and Section 4 contains the conclusions. 



2 Theoretical framework and calculation technique 



The construction of the effective lagrangian for the electroweak interactions 
satisfying the requirement of SU(2) x U(l) gauge invariance may be found in 
the literature [8]. It will be useful to mention here that in this approach the 
lagrangian can be written in the form 

C = Cym + C F + C S , (1) 



where Cym is the Yang-Mills lagrangian, 

£y« = -\w; u W^ - \b, v B^ , (2) 



Cf is the fermionic contribution and Cs is the scalar contribution. The stan- 
dard couplings between gauge bosons are given by the Yang-Mills lagrangian. 
The scalar sector is represented by means of the unitary matrix 

U = exp (i T -f) , (3) 



where v = 246 GeV and 7T; are the would-be Goldstone bosons. The gauge- 
invariant operators contributing to Cs can be constructed by taking the traces 
of the following building blocks: 



T 



w, 



■dnWv-dvWu + iglW^Wv] 



(4) 
(5) 
(6) 
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where D^U = dJJ + ^t^W^U — '^B^Ut 3 . At leading order, assuming the 
custodial symmetry, the following contribution to C s is found: 

C = V \r ((D^D^U) , (7) 



that gives the masses to the W ± and Z bosons. 

At next-to-leading order, that corresponds to dimension four, new couplings 
appear, and the independent gauge-invariant structures are multiplied by co- 
efficients that are model-dependent. At this order, trilinear and quadrilinear 
vertices are present. The operators that give rise to quadrilinear and not to tri- 
linear vertices, with the condition of CP conservation, are (using the standard 
notation): 



C 4 = a 4 (Tr(V,V u )) 2 

C 5 = a 5 (Tr(V^V»)) 2 

£ 6 = a 6 Tr(y^V v )Tr(TV»)Tr(TV») 

C 7 = a 7 Tr(y^)(Tr(TV u )) 2 

C w = a w l -{Tr{TV,)Tr{TV v )f . (8) 

The custodial symmetry is respected only by the operators £4 and £ 5 , and is 
violated by the other three operators, due to the presence of T. In the unitary 
gauge the above terms are given by: 



(9) 



+ ±.w+z*w;z» + -^(z.zn 2 ) 
£5 = « 5 ( (w+w-») 2 + ^w^w-^z v z v 

+ ik- iz " z " )2 ) (10) 

£ 6 = a„g 4 f-i-»; + z"»'-z" + -L(z„z'A (11) 

\ c w zc w / 

£7 = a 7 g 4 (±W+W-»Z„Z» + ^(Z^) 2 ] (12) 
\ c w Ac w ) 

C 10 = a 107 ^(Z^) 2 . (13) 
Zc w 
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As can be seen, these operators give the anomalous contributions to the quadri- 
linear gauge couplings involving the W ± and Z vector bosons, to be added 
to the standard ones that are contained in the Yang-Mills lagrangian. In par- 
ticular, it is easy to find, in the above formulas, anomalous contributions to 
the vertices AW (£ 4 and £ 5 ), WWZZ (£ 4 , £5, A and £ 7 ) and AZ (all the 
operators). 

In this paper, the numerical results of a phenomenological study on 6/ pro- 
cesses involving these anomalous couplings will be presented. These results 
have been obtained by means of a computer code already employed in other 
6/ analyses [9,10]. In this code the scattering amplitudes are calculated by the 
automatic algorithm ALPHA [11], and the Monte Carlo integration procedure 
is the result of an adaptation of the four-fermion codes HIGGSPV [12] and 
WWGENPV [13] to the 6/ calculations. For the present study, the lagrangian 
in eq. (1) with C s = C + £ 4 + £ 5 + C e + C 7 + C 10 has been implemented. 
The case in which the coefficients ctj are all equal to zero is considered as a 
reference model, that represents the limit of infinite Higgs mass in the SM 
at tree-level, and the deviations from such a model when the coefficients are 
non-vanishing have been studied. 

The input parameters are G M , M w and M z . The widths of the W and Z° 
bosons and all the couplings are calculated at tree-level. All the fermions are 
massless. For the propagators of the gauge bosons, the "fixed-width" scheme 
has been adopted. The reliability of this approach in 6/ calculations at energies 
of the order of the TeV has been discussed in detail in ref. [9]. 

The aim of the present work is to provide an analysis, that has never been 
done till now, of the impact of a complete electroweak 6/ calculation on the 
study of anomalous quartic couplings. The numerical results contained in the 
following section are a first illustrative application of the new implementa- 
tion of anomalous couplings in the code mentioned above, and the processes 
considered, although not giving a complete picture of the anomalous quartic 
coupling phenomenology, are relevant due to their sizeable cross-sections, as 
can be seen below. 



3 Physical processes and results 

The processes e + e~ — > 2q + 2q'v e & e , with q = u,c and q' = d,s, have been 
considered. This choice is motivated by the two objectives of having contri- 
butions from all possible quartic vertices involved in the anomalous terms 
under consideration, and of having no more than two neutrinos in the final 
state. The first point can be explained by considering the signature uuddu e u e : 
as can be seen in fig. 1, that shows the diagrams with a quadrilinear vertex, 
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Fig. 1. Diagrams with quadrilinear massive gauge boson vertices contributing to 
the process e + e~ — > uuddv e u e . Notice that the diagram on the right in the first row 
(with a AZ vertex) can be only anomalous, while the others have both standard and 
anomalous contributions. 

this signature includes contributions from the 4W, WW ZZ and 4Z vertices. 
From the point of view of the calculation, in order to obtain the correct sum 
over colours, with the version of ALPHA used here, that does not contain the 
colour degrees of freedom, it is necessary to apply a procedure analogous to 
the one already illustrated in ref. [10]. This procedure consists of the combi- 
nation, with proper weights, of the results from the different processes and 
is valid in the approximation of massless quarks and unit CKM matrix. The 
signatures udcsv e v e , uuddu e u e and uussu e u e must be combined according to 
the following formula: 



where the results directly provided by ALPHA are denoted by a, while the 
cross-sections including the colour degrees of freedom are denoted by a. The 
interested reader is referred to [10] for more details. 

The signatures csudv e V e , ccssv e v e and ccddv e v e are obtained from the previous 
ones by means of the simultaneous exchanges u <-> c and d <-> s, and, given 
the set of parameters adopted here, they give exactly the same result, so that 
they can be easily taken into account by including a factor of 2 in eq. (14). 

In principle, assuming that the 6-tagging technique is applied, and thus the 
b quark can be identified, a realistic calculation should take into account all 
the possible combinations of the remaining flavours, u, d, c and s that cannot 
be distinguished. These combinations include, in addition to those mentioned 



@uudd @udcs 



uuss 



+ (2N C -1)( 

^udcs "uuss )) , (14) 
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above, the final states uuccv e v e , uuuuv e v e ddssv e v e ddddv e v e and those ob- 
tained from these with the exchanges u <-> c and d «-> s. Moreover a sum over 
the neutrino flavours should be made. Such signatures have been neglected in 
this first analysis. Nevertheless, it will be possible to examine the most impor- 
tant qualitative features of the phenomenology of anomalous gauge couplings 
in Qf processes, while referring to further developments for a better accuracy 
from a quantitative point of view. 

In the following, the results of complete electroweak calculations at tree level 
are presented, where the above specified class of processes is considered. 

As can be seen in the scattering of real gauge bosons, as a consequence of the 
absence of the Higgs boson and of the gauge cancellations that ensure uni- 
tarity, the cross-sections at high energy can be expected to increase rapidly 
and to violate the unitarity bounds. The energy at which unitarity is violated 
depends on the values of the parameters c^: these should then be taken within 
some bounds if the unitarity condition has to be respected at a given energy. 
However, it is useful to make a first rough analysis of the dependence of the 
cross-sections on the anomalous parameters in a wide range, while the sensi- 
tivity to the anomalous couplings in a smaller range will be studied in a second 
step. Thus, in the first set of results, the cross-section has been evaluated at 
the two cm. energies of 500 GeV and 1 TeV, and allowing each of the coeffi- 
cients ccj to vary in the range (—0.2, 0.2), while keeping all the others equal to 
zero. A simple set of kinematical constraints has been adopted, by requiring 
the invariant masses of the "up-anti-up" and "down- ant i- down" quark pairs 
to be greater than 70 GeV (where "up" stands for u and c and "down" for d 
and s), so as to eliminate the soft photon contributions to pair production. 

In the plots of fig. 2 the cross-section is shown as a function of the various 
parameters ctj at the two energies of 500 GeV and 1 TeV. First of all, it 
should be observed that the effects of the anomalous couplings are greater at 
1 TeV than at 500 GeV, as expected, according to the above arguments on 
unitarity violation. As a consequence, the cross-section turns out to have a 
very strong growth with energy for large values of the anomalous coefficients. 
The growth of the cross-section with energy is in agreement with results for 
the scattering of real longitudinal gauge bosons obtained by means of the 
equivalence theorem [7]. Moreover it can be observed that the greatest effects 
are given by the coefficients and a$, while the coefficient a±o gives the 
smallest effect. To understand this fact it is useful to examine the expressions 
of the operators in eqs. (9-13), where it can be seen that £ 4 and £ 5 are the only 
terms where the vertex 4W appears. This vertex gives a greater enhancement 
with respect to the others due to the fact that the couplings of the W boson to 
the fermions are stronger than those of the Z boson. For the same reason, the 
operator C w has the smallest effect, since it contributes only to the 4Z vertex. 
Moreover the AZ vertex, as can be seen in fig. 1, can appear only in s-channel 
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Fig. 2. Cross-section in the Born approximation as a function of the anoma- 
lous coupling parameters. The invariant masses of the pairs of up-anti-up and 
down-anti-down quarks are greater than 70 GeV. 

diagrams, and has not the advantage of the t-channel growth at high energy. 

The parameters a^, a-j and a±o, that induce violation of the custodial symme- 
try, are more constrained by the radiative corrections to the p parameter with 
respect to and ot$. For this reason, the latter are more interesting and the 
remaining part of the study will be restricted to them. 

The behaviour of the cross-section as a function of the cm. energy is considered 
in fig. 3, where the energy range between 500 GeV and 1 TeV is studied for 
the case without anomalous couplings and for two choices of the anomalous 
parameters taken as examples, a 4 = 0.05 and a 5 = 0.05. The invariant masses 
of all the possible quark pairs are required to be above 60 GeV. The solid 
curves include initial-state radiation [14] (ISR) and beamstrahlung [15] (BS), 
while the dashed ones are in the Born approximation. The growth with energy 
is already present when the parameters ctj are equal to zero, and is enhanced 
when they are different from zero. The effect of ISR and BS is to reduce the 
effective cm. energy, and this explains the observed lowering of 15 — 20% in 
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Fig. 3. Integrated cross-section with initial-state radiation and beamstrahlung (solid 
lines) and in the Born approximation (dashed lines) as a function of the cm. energy. 
The result with all the coefficients ctj set to zero is compared with the cases 014 = 0.05 
(first plot) and 05 = 0.05 (second plot). 

the cross-sections with respect to the Born approximation. 

In order to study the sensitivity of some differential distributions to the anoma- 
lous couplings, two samples of events have then been generated, one with all 
the ctj equal to zero and the other with 0:4 = 0.1. The cm. energy is 1 TeV 
and the numbers of events correspond to a luminosity of 1000 fb _1 . As in 
fig. 3, all the quark pairs have an invariant mass greater than 60 GeV. A 
relatively large value of 0:4 has been chosen, so as to have sizeable effects: 
indeed, the objective of this analysis is to obtain clear indications on the cuts 
to be applied for enhancing the sensitivity to the anomalous couplings under 
study. Among the variables considered, the ones that turn out to be most 
sensitive to the parameter a 4 are denned in the following way. The invariant 
mass of the system of four jets is indicated as M(WW). This variable does 
not require any identification procedure. The other variables that have been 
studied are based instead on simple identification algorithms, to take into ac- 
count the fact that the quark flavours are not distinguishable. The W bosons 
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are reconstructed as the pairs of quarks q^q^ and quqi such that the quantity 
\m(q,iqj) 2 — + \m(qkqi) 2 — is minimized. The angle of one recon- 

structed W boson with respect to the beam axis is indicated as 6{W). For 
each event, the four jets are ordered in transverse momentum and are labelled 
with ji, . . . , j 4 , where j\ is the one with greatest p t . The invariant masses of 
pairs of jets are then considered, and one in particular, M(j^j4), that is the 
invariant mass of the pair of jets with lowest pt, is found to be most sensitive 
to the anomalous couplings under study. In fig. 4 the variables defined above 
are shown for the sample with a 4 = 0.1 (solid histograms) in comparison with 
the sample where all the anomalous parameters are set to zero (dashed his- 
tograms). In these plots, where the different total numbers of events for the 
two samples, due to the differences of the cross-sections, must be taken into 
account, significant deviations are found in the shapes of the distributions. 
This is due to the Lorentz structures of the anomalous vertices involved, that 
tend to populate the region of phase-space with high transverse momenta for 
the gauge bosons. In view of a realistic simulation, these distributions include 
the effects of ISR and BS. It has been verified that these effects do not intro- 
duce significant variations with respect to the Born approximation, as can be 
expected, since the variables considered do not involve the momenta of the 
neutrinos. 

On the contrary, the distribution of the missing mass, defined as M miss = 



V Pmissi where P m iss is the missing momentum, is strongly affected by ISR 
and BS. This variable, that in the Born approximation is the invariant mass 
of the neutrino pair, is shown in fig. 5: in the upper plot the case with ctj = 
is considered, while the lower plot refers to a 4 = 0.1. The solid histograms 
include ISR and BS, the dashed histograms are in the Born approximation. 
In both plots the peaks corresponding to the Z mass are almost cancelled by 
the effects of ISR and BS. 

Results very similar to those illustrated above are obtained for the coupling 
«5. On the basis of this analysis of distributions, suitable cuts have then 
been determined to enhance the sensitivity to the anomalous couplings a 4 
and a§. The results are given in fig. 6, where the cross-section in the Born 
approximation is shown as a function of the parameters 0:4 (upper plot) and 
a 5 (lower plot) in the range (—0.01,0.01). 

The cuts are defined as follows: M(WW) is greater than 420 GeV, the in- 
variant mass of the two jets with the lowest transverse momenta, m(j 3 j4), is 
greater than 80 GeV and the angle of one reconstructed W boson with re- 
spect to the beam axis satisfies cos 9(W) < 0.7. Moreover, the invariant mass 
of all the quark pairs is greater than 60 GeV. It has been found that neither 
the addition of a cut on the variable pt{W) nor its substitution to one of the 
above cuts can improve the results obtained. In the same plots, the la limits 
around the value of the cross-section for ctj = 0, with an experimental error 
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Vs=1 000 GeV 




Fig. 4. Comparison between two samples of events with 04 = 0.1 (solid histograms) 
and a>i = (dashed histograms) in the presence of initial-state-radiation and beam- 
strahlung. The samples are normalized to a luminosity of 1000 fb _1 , and have dif- 
ferent numbers of events, due to the different cross-sections . (a): invariant mass 
distribution of the system of four jets; (b): distribution of the transverse momen- 
tum of the reconstructed W boson; (c): angular distribution of the reconstructed 
W boson with respect to the beam axis; (d): invariant mass distribution of the pair 
of quarks with lowest transverse momenta. 

corresponding to a luminosity of 1000 fir 1 , are shown. It can be seen that the 
sensitivity to the parameters 0:4 and 0:5 turns out to be in the range between 
10 -3 and 10~ 2 . These conclusions are not modified by the introduction of ISR 
and BS, since the variables used in the cuts are not sensitive to these effects. 



4 Conclusions 



The analysis of anomalous dimension- four quartic gauge couplings involves 
processes with at least six fermions in the final state. A set of processes of 
this kind has been considered in this work. By using a Monte Carlo event 
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Fig. 5. Missing mass distribution. Upper plot: on = 0, lower plot: 04 = 0.1. Solid: 
ISR+BS; dashed: Born approximation. 

generator that has already been employed for other phenomenological studies 
on 6f physics, complete tree-level results have been obtained, in the context 
of the electroweak chiral lagrangian formalism. The five dimension-four oper- 
ators giving genuine quartic anomalous couplings have been considered, and a 
special attention has been devoted to the two custodial-symmetry conserving 
ones, usually indicated by £4 and £5 and involving the parameters 0:4 and a§. 
After studying the energy dependence of the cross-section and the effects of 
initial-state-radiation and beamstrahlung, a set of kinematical cuts has been 
considered in order to enhance the signals of the anomalous couplings. A sen- 
sitivity to the parameters «4 and a 5 in the range between 10~ 3 and 10~ 2 has 
been found at 1 TeV. 

The results discussed above can be seen as an example of the possibilities of 
study of these phenomena through a complete tree-level simulation of 6f pro- 
cesses. In order to make the calculations quantitatively more accurate, some 
simplifications that have been adopted here should be eliminated: namely the 
class of signatures to be taken into account should include a full sum over the 
indistinguishable neutrino and quark flavour combinations in the final state, 
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Vs = 1 TeV 
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Fig. 6. Integrated cross-section in the Born approximation at a cm. energy of 1 TeV, 
as a function of 04 (first plot) and 05 (second plot) with a set of kinematical cuts 
studied to enhance the sensitivity to the anomalous couplings. The horizontal solid 
lines are the la bounds around the a\ = value (dashed line), with the experimental 
uncertainty corresponding to a luminosity of 1000 fb _1 . 



and the effect of QCD backgrounds should be considered. In particular the lat- 
ter improvement can be achieved by means of the last version of ALPHA [16], 
that includes the QCD lagrangian. Moreover, the role of polarization of the 
initial electrons and positrons could be investigated. These points will be the 
the subject of future developments. 
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